Exploring the pediatric nasopharyngeal bacterial microbiota with culture-based MALDI-TOF mass spectrometry and targeted metagenomic sequencing.

The nasopharynx is an important reservoir of disease-associated and antimicrobial-resistant bacterial species. This proof-of-concept study assessed the utility of a combined culture, matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS), and targeted metagenomic sequencing workflow for the study of the pediatric nasopharyngeal bacterial microbiota. Nasopharyngeal swabs and clinical metadata were collected from Cambodian children during a hospital outpatient visit and then biweekly for 12 weeks. Swabs were cultured on chocolate and blood-gentamicin agar, and all colony morphotypes were identified by MALDI-TOF MS. Metagenomic sequencing was done on a scrape of all colonies from a chocolate agar culture and processed using the mSWEEP pipeline. One hundred one children were enrolled, yielding 620 swabs. MALDI-TOF MS identified 106 bacterial species/40 genera: 20 species accounted for 88.5% (2,190/2,474) of isolates. Colonization by Moraxella catarrhalis (92.1% of children on ≥1 swab), Haemophilus influenzae (87.1%), and Streptococcus pneumoniae (83.2%) was particularly common. In S. pneumoniae-colonized children, a median of two serotypes [inter-quartile range (IQR) 1-2, range 1-4] was detected. For the 21 bacterial species included in the mSWEEP database and identifiable by MALDI-TOF, detection by culture + MALDI-TOF MS and culture + mSWEEP was highly concordant with a median species-level agreement of 96.9% (IQR 86.8%-98.8%). mSWEEP revealed highly dynamic lineage-level colonization patterns for S. pneumoniae which were quite different to those for S. aureus. A combined culture, MALDI-TOF MS, targeted metagenomic sequencing approach for the exploration of the young child nasopharyngeal microbiome was technically feasible, and each component yielded complementary data. IMPORTANCE: The human upper respiratory tract is an important source of disease-causing and antibiotic-resistant bacteria. However, understanding the interactions and stability of these bacterial populations is technically challenging. We used a combination of approaches to determine colonization patterns over a 3-month period in 101 Cambodian children. The combined approach was feasible to implement, and each component gave complementary data to enable a better understanding of the complex patterns of bacterial colonization.

Epistasis, core-genome disharmony, and adaptation in recombining bacteria.

Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome's co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background-facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche.IMPORTANCEGenetic exchange among bacteria shapes the microbial world. From the acquisition of antimicrobial resistance genes to fundamental questions about the nature of bacterial species, this powerful evolutionary force has preoccupied scientists for decades. However, the mixing of genes between species rests on a paradox: 0n one hand, promoting adaptation by conferring novel functionality; on the other, potentially introducing disharmonious gene combinations (negative epistasis) that will be selected against. Taking an interdisciplinary approach to analyze natural populations of the enteric bacteria Campylobacter, an ideal example of long-range admixture, we demonstrate that genes can independently transfer across species boundaries and rejoin in functional networks in a recipient genome. The positive impact of two-gene interactions appears to be adaptive by expanding metabolic capacity and facilitating niche shifts through interspecific hybridization. This challenges conventional ideas and highlights the possibility of multiple-step evolution of multi-gene traits by interspecific introgression.

Estimating Between Country Migration in Pneumococcal Populations.

Streptococcus pneumoniae (the pneumococcus) is a globally distributed, human obligate opportunistic bacterial pathogen which, although often carried commensally, is also a significant cause of invasive disease. Apart from multi-drug resistant and virulent clones, the rate and direction of pneumococcal dissemination between different countries remains largely unknown. The ability for the pneumococcus to take a foothold in a country depends on existing population configuration, the extent of vaccine implementation, as well as human mobility since it is a human obligate bacterium. To shed light on its international movement, we used extensive genome data from the Global Pneumococcal Sequencing (GPS) project and estimated migration parameters between multiple countries in Africa. Data on allele frequencies of polymorphisms at housekeeping-like loci for multiple different lineages circulating in the populations of South Africa, Malawi, Kenya, and The Gambia were used to calculate the fixation index (Fst) between countries. We then further used these summaries to fit migration coalescent models with the likelihood-free inference algorithms available in the ELFI software package. Synthetic data were additionally used to validate the inference approach. Our results demonstrate country-pair specific migration patterns and heterogeneity in the extent of migration between different lineages. Our approach demonstrates that coalescent models can be effectively used for inferring migration rates for bacterial species and lineages provided sufficiently granular population genomics surveillance data. Further it can demonstrate the connectivity of respiratory disease agents between countries to inform intervention policy in the longer term.

PlasmidEC and gplas2: an optimized short-read approach to predict and reconstruct antibiotic resistance plasmids in Escherichia coli.

Accurate reconstruction of Escherichia coli antibiotic resistance gene (ARG) plasmids from Illumina sequencing data has proven to be a challenge with current bioinformatic tools. In this work, we present an improved method to reconstruct E. coli plasmids using short reads. We developed plasmidEC, an ensemble classifier that identifies plasmid-derived contigs by combining the output of three different binary classification tools. We showed that plasmidEC is especially suited to classify contigs derived from ARG plasmids with a high recall of 0.941. Additionally, we optimized gplas, a graph-based tool that bins plasmid-predicted contigs into distinct plasmid predictions. Gplas2 is more effective at recovering plasmids with large sequencing coverage variations and can be combined with the output of any binary classifier. The combination of plasmidEC with gplas2 showed a high completeness (median=0.818) and F1-Score (median=0.812) when reconstructing ARG plasmids and exceeded the binning capacity of the reference-based method MOB-suite. In the absence of long-read data, our method offers an excellent alternative to reconstruct ARG plasmids in E. coli.

Co-culturing with Streptococcus anginosus alters Staphylococcus aureus transcriptome when exposed to tonsillar cells.

Introduction: Improved understanding of Staphylococcus aureus throat colonization in the presence of other co-existing microbes is important for mapping S. aureus adaptation to the human throat, and recurrence of infection. Here, we explore the responses triggered by the encounter between two common throat bacteria, S. aureus and Streptococcus anginosus, to identify genes in S. aureus that are important for colonization in the presence of human tonsillar epithelial cells and S. anginosus, and further compare this transcriptome with the genes expressed in S. aureus as only bacterium. Methods: We performed an in vitro co-culture experiment followed by RNA sequencing to identify interaction-induced transcriptional alterations and differentially expressed genes (DEGs), followed by gene enrichment analysis. Results and discussion: A total of 332 and 279 significantly differentially expressed genes with p-value < 0.05 and log2 FoldChange (log2FC) ≥ |2| were identified in S. aureus after 1 h and 3 h co-culturing, respectively. Alterations in expression of various S. aureus survival factors were observed when co-cultured with S. anginosus and tonsillar cells. The serine-aspartate repeat-containing protein D (sdrD) involved in adhesion, was for example highly upregulated in S. aureus during co-culturing with S. anginosus compared to S. aureus grown in the absence of S. anginosus, especially at 3 h. Several virulence genes encoding secreted proteins were also highly upregulated only when S. aureus was co-cultured with S. anginosus and tonsillar cells, and iron does not appear to be a limiting factor in this environment. These findings may be useful for the development of interventions against S. aureus throat colonization and could be further investigated to decipher the roles of the identified genes in the host immune response in context of a throat commensal landscape.

Modulation of multidrug-resistant clone success in Escherichia coli populations: a longitudinal, multi-country, genomic and antibiotic usage cohort study.

BACKGROUND: The effect of antibiotic usage on the success of multidrug-resistant (MDR) clones in a population remains unclear. With this genomics-based molecular epidemiology study, we aimed to investigate the contribution of antibiotic use to Escherichia coli clone success, relative to intra-strain competition for colonisation and infection. METHODS: We sequenced all the available E coli bloodstream infection isolates provided by the British Society for Antimicrobial Chemotherapy (BSAC) from 2012 to 2017 (n=718) and combined these with published data from the UK (2001-11; n=1090) and Norway (2002-17; n=3254). Defined daily dose (DDD) data from the European Centre for Disease Prevention and Control (retrieved on Sept 21, 2021) for major antibiotic classes (ß-lactam, tetracycline, macrolide, sulfonamide, quinolone, and non-penicillin ß-lactam) were used together with sequence typing, resistance profiling, regression analysis, and non-neutral Wright-Fisher simulation-based modelling to enable systematic comparison of resistance levels, clone success, and antibiotic usage between the UK and Norway. FINDINGS: Sequence type (ST)73, ST131, ST95, and ST69 accounted for 892 (49·3%) of 1808 isolates in the BSAC collection. In the UK, the proportion of ST69 increased between 2001-10 and 2011-17 (p=0·0004), whereas the proportions of ST73 and ST95 did not vary between periods. ST131 expanded quickly after its emergence in 2003 and its prevalence remained consistent throughout the study period (apart from a brief decrease in 2009-10). The extended-spectrum ß-lactamase (ESBL)-carrying, globally disseminated MDR clone ST131-C2 showed overall greater success in the UK (154 [56·8%] of 271 isolates in 2003-17) compared with Norway (51 [18·3%] of 278 isolates in 2002-17; p<0·0001). DDD data indicated higher total use of antimicrobials in the UK, driven mainly by the class of non-penicillin ß-lactams, which were used between 2·7-times and 5·1-times more in the UK per annum (ratio mean 3·7 [SD 0·8]). This difference was associated with the higher success of the MDR clone ST131-C2 (pseudo-R2 69·1%). A non-neutral Wright-Fisher model replicated the observed expansion of non-MDR and MDR sequence types under higher DDD regimes. INTERPRETATION: Our study indicates that resistance profiles of contemporaneously successful clones can vary substantially, warranting caution in the interpretation of correlations between aggregate measures of resistance and antibiotic usage. Our study further suggests that in countries with low-to-moderate use of antibiotics, such as the UK and Norway, the extent of non-penicillin ß-lactam use modulates rather than determines the success of widely disseminated MDR ESBL-carrying E coli clones. Detailed understanding of underlying causal drivers of success is important for improved control of resistant pathogens. FUNDING: Trond Mohn Foundation, Marie Sklodowska-Curie Actions, European Research Council, Royal Society, and Wellcome Trust.

BurkHostGEN: a study protocol for evaluating variations in the Burkholderia pseudomallei and host genomes associated with melioidosis infection.

Background: Melioidosis is a frequently fatal disease caused by an environmental bacterium Burkholderia pseudomallei. The disease is prevalent in northeast Thailand, particularly among rice field farmers who are at risk of bacterial exposure through contact with contaminated soil and water. However, not all exposure results in disease, and infection can manifest diverse outcomes. We postulate that genetic factors, whether from the bacterium, the host or the combination of both, may influence disease outcomes. To address this hypothesis, we aim to collect, sequence, and analyse genetic data from melioidosis patients and controls, along with isolates of B. pseudomallei obtained from patients. Additionally, we will study the metagenomics of the household water supply for both patients and controls, including the presence of B. pseudomallei. Methods: BurkHostGEN is an ongoing observational study being conducted at Sunpasitthiprasong Hospital, Ubon Ratchathani, Thailand. We are obtaining consent from 600 melioidosis patients and 700 controls, spanning both sexes, to collect 1 mL of blood for host DNA analysis, 3 mL of blood for RNA analysis, as well as 5 L of household water supply for metagenomic analysis. Additionally, we are isolating B. pseudomallei from the melioidosis patients to obtain bacterial DNA. This comprehensive approach will allow us to identify B. pseudomallei and their paired host genetic factors associated with disease acquisition and severity. Ethical approvals have been obtained for BurkHostGEN. Host and bacterial genetic data will be uploaded to European Genome-Phenome Archive (EGA) and European Nucleotide Archive (ENA), respectively. Conclusions: BurkHostGEN holds the potential to discover bacterial and host genetic factors associated with melioidosis infection and severity of illness. It can also support various study designs, including biomarker validation, disease pathogenesis, and epidemiological analysis not only for melioidosis but also for other infectious diseases.

Multidrug-resistant E. coli encoding high genetic diversity in carbohydrate metabolism genes displace commensal E. coli from the intestinal tract.

Extra-intestinal pathogenic Escherichia coli (ExPEC) can cause a variety of infections outside of the intestine and are a major causative agent of urinary tract infections. Treatment of these infections is increasingly frustrated by antimicrobial resistance (AMR) diminishing the number of effective therapies available to clinicians. Incidence of multidrug resistance (MDR) is not uniform across the phylogenetic spectrum of E. coli. Instead, AMR is concentrated in select lineages, such as ST131, which are MDR pandemic clones that have spread AMR globally. Using a gnotobiotic mouse model, we demonstrate that an MDR E. coli ST131 is capable of out-competing and displacing non-MDR E. coli from the gut in vivo. This is achieved in the absence of antibiotic treatment mediating a selective advantage. In mice colonised with non-MDR E. coli strains, challenge with MDR E. coli either by oral gavage or co-housing with MDR E. coli colonised mice results in displacement and dominant intestinal colonisation by MDR E. coli ST131. To investigate the genetic basis of this superior gut colonisation ability by MDR E. coli, we assayed the metabolic capabilities of our strains using a Biolog phenotypic microarray revealing altered carbon metabolism. Functional pangenomic analysis of 19,571 E. coli genomes revealed that carriage of AMR genes is associated with increased diversity in carbohydrate metabolism genes. The data presented here demonstrate that independent of antibiotic selective pressures, MDR E. coli display a competitive advantage to colonise the mammalian gut and points to a vital role of metabolism in the evolution and success of MDR lineages of E. coli via carriage and spread.

The population-level impact of Enterococcus faecalis genetics on intestinal colonization and extraintestinal infection.

IMPORTANCE: Enterococcus faecalis causes life-threatening invasive hospital- and community-associated infections that are usually associated with multidrug resistance globally. Although E. faecalis infections cause opportunistic infections typically associated with antibiotic use, immunocompromised immune status, and other factors, they also possess an arsenal of virulence factors crucial for their pathogenicity. Despite this, the relative contribution of these virulence factors and other genetic changes to the pathogenicity of E. faecalis strains remain poorly understood. Here, we investigated whether specific genomic changes in the genome of E. faecalis isolates influence its pathogenicity-infection of hospitalized and nonhospitalized individuals and the propensity to cause extraintestinal infection and intestinal colonization. Our findings indicate that E. faecalis genetics partially influence the infection of hospitalized and nonhospitalized individuals and the propensity to cause extraintestinal infection, possibly due to gut-to-bloodstream translocation, highlighting the potential substantial role of host and environmental factors, including gut microbiota, on the opportunistic pathogenic lifestyle of this bacterium.

Population structure and antimicrobial resistance among Klebsiella isolates sampled from human, animal, and environmental sources in Ghana: a cross-sectional genomic One Health study.

BACKGROUND: One Health approaches to address the increasing threat of antimicrobial resistance (AMR) are gaining attention. However, data on the distribution and movement of bacteria and their AMR-associated genes between clinical and non-clinical sources are scarce, especially from low-income and middle-income countries. We aimed to analyse Klebsiella isolates from various sources in Ghana and compare the prevalence of AMR with datasets from two other countries. METHODS: We conducted a cross-sectional genomic One Health study. Multiple clinical, environmental, and animal sources were sampled from 78 locations (eg, hospitals, residential areas, and farms) in and around Tamale, Ghana. Clinical samples were collected through routine screening and in cases of suspected infection between March 15 and Sept 15, 2019, and samples from the wider environment were collected during a dedicated sampling effort between the dates of Aug 19, 2018, and Sept 26, 2019. Sampling locations were approximately evenly distributed from the centre of the city and steadily outwards to capture both rural and urban locations. Samples with positive growth for Klebsiella were included. Isolates of Klebsiella were obtained from the samples using Simmons citrate agar medium and characterised by antimicrobial susceptibility testing and whole-genome sequencing. A comparative analysis with Klebsiella population surveys from Pavia, Italy, and Tromsø, Norway, was performed. AMR-associated and virulence genes were detected, and the population distribution of these genes was studied. FINDINGS: Of 957 samples collected around Tamale, Ghana, 620 were positive for Klebsiella spp. 573 Klebsiella isolates were successfully sequenced, of which 370 were Klebsiella pneumoniae. Only two hospital isolates were carbapenem-resistant. Extended-spectrum ß-lactamase (ESBL) genes were relatively common among the Ghanaian clinical isolates but rare in the environmental samples. Prevalence of ESBL genes in human-hospital disease samples was 64% (14 of 22 isolates) in Ghana and 44% (four of nine isolates) in Italy, and prevalence in human-hospital carriage samples was 7% (eight of 107) in Ghana and 13% (54 of 428) in Italy; the prevalence was higher in human-hospital disease samples than in human-hospital carriage samples in both countries, and prevalence across both samples in both countries was higher than in Norway. Ghanaian isolates showed evidence of high recombination rates (recombination events compared with point mutations [r/m] 9·455) and a considerable accessory gene overlap with isolates from Italy and Norway. INTERPRETATION: Although several AMR-associated gene classes were observed relatively frequently in non-clinical sources, ESBL, carbapenemase, and virulence genes were predominantly present only in hospital samples. These results suggest that interventions should be focused on clinical settings to have the greatest effect on the prevalence and dissemination of AMR-associated genes. FUNDING: European Research Council (742158), Academy of Finland EuroHPC grant, Trond Mohn Foundation (BATTALION grant), and Wellcome Trust.

Unraveling Enterococcus susceptibility to quaternary ammonium compounds: genes, phenotypes, and the impact of environmental conditions.

Quaternary ammonium compounds (QACs) have been extensively used in the community, healthcare facilities, and food chain, in concentrations between 20 and 30,000 mg/L. Enterococcus faecalis and Enterococcus faecium are ubiquitous in these settings and are recognized as nosocomial pathogens worldwide, but QACs' activity against strains from diverse epidemiological and genomic backgrounds remained largely unexplored. We evaluated the role of Enterococcus isolates from different sources, years, and clonal lineages as hosts of QACs tolerance genes and their susceptibility to QACs in optimal, single-stress and cross-stress growth conditions. Only 1% of the Enterococcus isolates included in this study and 0.5% of publicly available Enterococcus genomes carried qacA/B, qacC, qacG, qacJ, qacZ, qrg, bcrABC or oqxAB genes, shared with >60 species of Bacillota, Pseudomonadota, Actinomycetota, or Spirochaetota. These genes were generally found within close proximity of antibiotics and/or metals resistance genes. The minimum inhibitory concentrations (MIC) and minimum bactericidal concentrations (MBC) of benzalkonium chloride (BC) and didecyldimethylammonium chloride ranged between 0.5 and 4 mg/L (microdilution: 37°C/20 h/pH = 7/aerobiosis) for 210 E. faecalis and E. faecium isolates (two isolates carrying qacZ). Modified growth conditions (e.g., 22°C/pH = 5) increased MICBC/MBCBC (maximum of eightfold and MBCBC = 16 mg/L) and changed bacterial growth kinetics under BC toward later stationary phases in both species, including in isolates without QACs tolerance genes. In conclusion, Enterococcus are susceptible to in-use QACs concentrations and rarely carry QACs tolerance genes. However, their potential gene exchange with different microbiota, the decreased susceptibility to QACs under specific environmental conditions, and the presence of subinhibitory QACs concentrations in various settings may contribute to the selection of particular strains and, thus, require a One Health strategy to maintain QACs effectiveness. IMPORTANCE Despite the increasing use of quaternary ammonium compounds (QACs), the susceptibility of pathogens to these antimicrobials remains largely unknown. Enterococcus faecium and Enterococcus faecalis are susceptible to in-use QACs concentrations and are not main hosts of QACs tolerance genes but participate in gene transfer pathways with diverse bacterial taxa exposed to these biocides. Moreover, QACs tolerance genes often share the same genetic contexts with antibiotics and/or metals resistance genes, raising concerns about potential co-selection events. E. faecium and E. faecalis showed increased tolerance to benzalkonium chloride under specific environmental conditions (22°C, pH = 5), suggesting that strains might be selected in settings where they occur along with subinhibitory QACs concentrations. Transcriptomic studies investigating the cellular mechanisms of Enterococcus adaptation to QACs tolerance, along with longitudinal metadata analysis of tolerant populations dynamics under the influence of diverse environmental factors, are essential and should be prioritized within a One Health strategy.

Mge-cluster: a reference-free approach for typing bacterial plasmids.

Extrachromosomal elements of bacterial cells such as plasmids are notorious for their importance in evolution and adaptation to changing ecology. However, high-resolution population-wide analysis of plasmids has only become accessible recently with the advent of scalable long-read sequencing technology. Current typing methods for the classification of plasmids remain limited in their scope which motivated us to develop a computationally efficient approach to simultaneously recognize novel types and classify plasmids into previously identified groups. Here, we introduce mge-cluster that can easily handle thousands of input sequences which are compressed using a unitig representation in a de Bruijn graph. Our approach offers a faster runtime than existing algorithms, with moderate memory usage, and enables an intuitive visualization, classification and clustering scheme that users can explore interactively within a single framework. Mge-cluster platform for plasmid analysis can be easily distributed and replicated, enabling a consistent labelling of plasmids across past, present, and future sequence collections. We underscore the advantages of our approach by analysing a population-wide plasmid data set obtained from the opportunistic pathogen Escherichia coli, studying the prevalence of the colistin resistance gene mcr-1.1 within the plasmid population, and describing an instance of resistance plasmid transmission within a hospital environment.

Evolutionary and functional history of the Escherichia coli K1 capsule.

Escherichia coli is a leading cause of invasive bacterial infections in humans. Capsule polysaccharide has an important role in bacterial pathogenesis, and the K1 capsule has been firmly established as one of the most potent capsule types in E. coli through its association with severe infections. However, little is known about its distribution, evolution and functions across the E. coli phylogeny, which is fundamental to elucidating its role in the expansion of successful lineages. Using systematic surveys of invasive E. coli isolates, we show that the K1-cps locus is present in a quarter of bloodstream infection isolates and has emerged in at least four different extraintestinal pathogenic E. coli (ExPEC) phylogroups independently in the last 500 years. Phenotypic assessment demonstrates that K1 capsule synthesis enhances E. coli survival in human serum independent of genetic background, and that therapeutic targeting of the K1 capsule re-sensitizes E. coli from distinct genetic backgrounds to human serum. Our study highlights that assessing the evolutionary and functional properties of bacterial virulence factors at population levels is important to better monitor and predict the emergence of virulent clones, and to also inform therapies and preventive medicine to effectively control bacterial infections whilst significantly lowering antibiotic usage.

The highly diverse plasmid population found in Escherichia coli colonizing travellers to Laos and its role in antimicrobial resistance gene carriage.

Increased colonization by antimicrobial-resistant organisms is closely associated with international travel. This study investigated the diversity of mobile genetic elements involved with antimicrobial resistance (AMR) gene carriage in extended-spectrum beta-lactamase (ESBL)-producing Escherichia coli that colonized travellers to Laos. Long-read sequencing was used to reconstruct complete plasmid sequences from 48 isolates obtained from the daily stool samples of 23 travellers over a 3 week period. This method revealed a collection of 105 distinct plasmids, 38.1 % (n=40) of which carried AMR genes. The plasmids in this population were diverse, mostly unreported and included 38 replicon types, with F-type plasmids (n=23) the most prevalent amongst those carrying AMR genes. Fine-scale analysis of all plasmids identified numerous AMR gene contexts and emphasized the importance of IS elements, specifically members of the IS6/IS26 family, in the evolution of complex multidrug resistance regions. We found a concerning convergence of ESBL and colistin resistance determinants, with three plasmids from two different F-type lineages carrying bla CTX-M and mcr genes. The extensive diversity seen here highlights the worrying probability that stable new vehicles for AMR will evolve in E. coli populations that can disseminate internationally through travel networks.

Challenges in prokaryote pangenomics.

Horizontal gene transfer (HGT) and the resulting patterns of gene gain and loss are a fundamental part of bacterial evolution. Investigating these patterns can help us to understand the role of selection in the evolution of bacterial pangenomes and how bacteria adapt to a new niche. Predicting the presence or absence of genes can be a highly error-prone process that can confound efforts to understand the dynamics of horizontal gene transfer. This review discusses both the challenges in accurately constructing a pangenome and the potential consequences errors can have on downstream analyses. We hope that by summarizing these issues researchers will be able to avoid potential pitfalls, leading to improved bacterial pangenome analyses.

Inferring the heritability of bacterial traits in the era of machine learning.

Quantification of heritability is a fundamental desideratum in genetics, which allows an assessment of the contribution of additive genetic variation to the variability of a trait of interest. The traditional computational approaches for assessing the heritability of a trait have been developed in the field of quantitative genetics. However, the rise of modern population genomics with large sample sizes has led to the development of several new machine learning-based approaches to inferring heritability. In this article, we systematically summarize recent advances in machine learning which can be used to infer heritability. We focus on an application of these methods to bacterial genomes, where heritability plays a key role in understanding phenotypes such as antibiotic resistance and virulence, which are particularly important due to the rising frequency of antimicrobial resistance. By designing a heritability model incorporating realistic patterns of genome-wide linkage disequilibrium for a frequently recombining bacterial pathogen, we test the performance of a wide spectrum of different inference methods, including also GCTA. In addition to the synthetic data benchmark, we present a comparison of the methods for antibiotic resistance traits for multiple bacterial pathogens. Insights from the benchmarking and real data analyses indicate a highly variable performance of the different methods and suggest that heritability inference would likely benefit from tailoring of the methods to the specific genetic architecture of the target organism. Availability and implementation: The R codes and data used in the numerical experiments are available at: https://github.com/tienmt/her_MLs.

The COVID-19 Pandemic Sparked Off a Large-Scale Outbreak of Carbapenem-Resistant Acinetobacter baumannii from the Endemic Strains at an Italian Hospital.

Acinetobacter baumannii is a nosocomial pathogen that poses a serious threat due to the rise of incidence of multidrug-resistant (MDR) strains. During the COVID-19 pandemic, MDR A. baumannii clones have caused several outbreaks worldwide. Here, we describe a detailed investigation of an MDR A. baumannii outbreak that occurred at Policlinico San Matteo (Pavia, Italy). A total of 96 A. baumannii strains, isolated between January and July 2020 from 41 inpatients (both SARS-CoV-2 positive and negative) in different wards, were characterized by phenotypic and genomic analyses combining Illumina and Nanopore sequencing. Antibiotic susceptibility testing revealed that all isolates were resistant to carbapenems, and the sequence analysis attributed this to the carbapenemase gene blaOXA-23. Virulence factor screening unveiled that all strains carried determinants for biofilm formation, while plasmid analysis revealed the presence of two plasmids, one of which was ~100 kbp long and encoded a phage sequence. A core genome-based phylogeny was inferred to integrate outbreak strain genomes with background genomes from public databases and the local surveillance program. All strains belonged to the globally disseminated sequence type 2 (ST2) clone and were mainly divided into two clades. Isolates from the outbreak clustered with surveillance isolates from 2019, suggesting that the outbreak was caused by two strains that were already circulating in the hospital before the start of the pandemic. The intensive spread of A. baumannii in the hospital was enhanced by the extreme emergency situation of the first COVID-19 pandemic wave that resulted in reduced attention to infection prevention and control practices. IMPORTANCE The COVID-19 pandemic, especially during the first wave, posed a great challenge to the hospital management and generally promoted nosocomial pathogen dissemination. MDR A. baumannii can easily spread and persist for a long time on surfaces, causing outbreaks in health care settings. Infection prevention and control practices, epidemiological surveillance, and microbiological screening are fundamental in order to control such outbreaks. Here, we sequenced the genomes of 96 isolates from an outbreak of MDR A. baumannii strains using both short- and long-read technology in order to reconstruct the outbreak events in fine detail. The sequence data demonstrated that two endemic clones of MDR A. baumannii were the source of this large hospital outbreak during the first COVID-19 pandemic wave, confirming the effect of COVID-19 emergency disrupting the protection provided by the use of the standard prevention procedures.

On predictive inference for intractable models via approximate Bayesian computation.

Approximate Bayesian computation (ABC) is commonly used for parameter estimation and model comparison for intractable simulator-based statistical models whose likelihood function cannot be evaluated. In this paper we instead investigate the feasibility of ABC as a generic approximate method for predictive inference, in particular, for computing the posterior predictive distribution of future observations or missing data of interest. We consider three complementary ABC approaches for this goal, each based on different assumptions regarding which predictive density of the intractable model can be sampled from. The case where only simulation from the joint density of the observed and future data given the model parameters can be used for inference is given particular attention and it is shown that the ideal summary statistic in this setting is minimal predictive sufficient instead of merely minimal sufficient (in the ordinary sense). An ABC prediction approach that takes advantage of a certain latent variable representation is also investigated. We additionally show how common ABC sampling algorithms can be used in the predictive settings considered. Our main results are first illustrated by using simple time-series models that facilitate analytical treatment, and later by using two common intractable dynamic models. Supplementary Information: The online version contains supplementary material available at 10.1007/s11222-022-10163-6.

Robust analysis of prokaryotic pangenome gene gain and loss rates with Panstripe.

Horizontal gene transfer (HGT) plays a critical role in the evolution and diversification of many microbial species. The resulting dynamics of gene gain and loss can have important implications for the development of antibiotic resistance and the design of vaccine and drug interventions. Methods for the analysis of gene presence/absence patterns typically do not account for errors introduced in the automated annotation and clustering of gene sequences. In particular, methods adapted from ecological studies, including the pangenome gene accumulation curve, can be misleading as they may reflect the underlying diversity in the temporal sampling of genomes rather than a difference in the dynamics of HGT. Here, we introduce Panstripe, a method based on generalized linear regression that is robust to population structure, sampling bias, and errors in the predicted presence/absence of genes. We show using simulations that Panstripe can effectively identify differences in the rate and number of genes involved in HGT events, and illustrate its capability by analyzing several diverse bacterial genome data sets representing major human pathogens.

Enterococcus spp. from chicken meat collected 20 years apart overcome multiple stresses occurring in the poultry production chain: Antibiotics, copper and acids.

Poultry meat has been a vehicle of antibiotic resistant bacteria and genes. Yet, the diversity of selective pressures associated with their maintenance in the poultry-production chain remains poorly explored. We evaluated the susceptibility of Enterococcus spp. from chicken meat collected 20 years apart to antibiotics, metals, acidic pH and peracetic acid-PAA. Contemporary chicken-meat samples (n = 53 batches, each including a pool of neck skin from 10 single carcasses) were collected in a slaughterhouse facility using PAA as disinfectant (March-August 2018, North of Portugal). Broilers were raised in intensive farms (n = 29) using CuSO4 and organic acids as feed additives. Data were compared with that of 67 samples recovered in the same region during 1999-2001. All 2018 samples had multidrug resistant-MDR isolates, with >45 % carrying Enterococcus faecalis, Enterococcus faecium or Enterococcus gallinarum resistant to tetracycline, erythromycin, ampicillin, quinupristin-dalfopristin, ciprofloxacin, chloramphenicol or aminoglycosides. Resistance rates were similar (P > 0.05) to those of 1999-2001 samples for all but five antibiotics. The decrease of samples carrying vancomycin-resistant isolates from 46 % to 0 % between 1999-2001 and 2018 was the most striking difference. Isolates from both periods were similarly susceptible to acid pH [minimum-growth pH (4.5-5.0), minimum-survival pH (3.0-4.0)] and to PAA (MIC90 = 100-120 mg/L/MBC90 = 140-160 mg/L; below concentrations used in slaughterhouse). Copper tolerance genes (tcrB and/or cueO) were respectively detected in 21 % and 4 % of 2018 and 1999-2001 samples. The tcrB gene was only detected in E. faecalis (MICCuSO4 > 12 mM), and their genomes were compared with other international ones of chicken origin (PATRIC database), revealing a polyclonal population and a plasmid or chromosomal location for tcrB. The tcrB plasmids shared diverse genetic modules, including multiple antimicrobial resistance genes (e.g. to tetracyclines, chloramphenicol, macrolide-lincosamide-streptogramin B-MLSB, aminoglycosides, bacitracin, coccidiostats). When in chromosome, the tcrB gene was co-located closely to merA (mercury) genes. Chicken meat remains an important vehicle of MDR Enterococcus spp. able to survive under diverse stresses (e.g. copper, acid) potentially contributing to these bacteria maintenance and flux among animal-environment-humans.